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Abstract 

Constraints on the core temperature (T c ) of the Sun and on neutrino-oscillation 
parameters are obtained from the existing solar neutrino data, including the recent 
GALLEX and Kamiokande III results. (1) A purely astrophysical solution to the solar 
neutrino problem is strongly disfavored by the data: the best fit in a cooler Sun model 
requires an 8% reduction in T c , but the x 2 test rejects this hypothesis at 99.99% C.L. 
This suggests neutrino physics (mass, mixings) that goes beyond the standard elec- 
troweak theory. (2) Assuming the Standard Solar Model (SSM) and matter-enhanced 
neutrino oscillations, the MSW parameters are constrained to two small regions: non- 
adiabatic oscillations with Am 2 = (0.3 - 1.2) x 1(T 5 eV 2 , sin 2 26 = (0.4 - 1.5) x 10" 2 , 
or large mixing-angle oscillations with Am 2 = (0.3 — 4) x 10 -5 eV 2 , sin 2 26 = 0.5 — 0.9. 
The non-adiabatic solution gives a considerably better fit. For v e oscillations into ster- 
ile neutrinos, the allowed region (90%) is constrained to non-adiabatic oscillations. As 
long as the SSM is assumed, the neutrino mixing angles are at least four times larger, 
or considerably smaller, than the corresponding quark mixing angles. (3) Allowing 
both MSW oscillations and a non-standard core temperature, a) the experiments de- 
termine the core temperature at the 5% level: T c = 1.03lg;o5 (90% C.L.) relative to 
the SSM, which is consistent with the SSM prediction, b) When T c is used as a free 
parameter, the allowed MSW region is broadened: a cooler Sun (T c = 0.95) allows 
Am 2 , sin 2 26 implied by the supersymmetric SO(10) grand unified theory (GUT), while 
a warmer Sun (T c = 1.05) extends the allowed parameter space into values suggested 
by intermediate-scale SO (10) GUTs, for which the v T may be cosmologically relevant. 
Superstring-inspired models are consistent with all solutions. (4) From the narrowed pa- 
rameter space, we predict the neutrino spectral shape which should be observed in the 
Sudbury Neutrino Observatory (SNO). Expected rates for SNO, Super-Kamiokande, 
and BOREXINO are also discussed. Throughout the calculation we use the Bahcall- 
Pinsonneault SSM (1992) with helium diffusion, and include nuclear and astrophysical 
uncertainties in a simplified, but physically transparent way. 

' Present address: Fermi National Accelerator Laboratory, P.O. Box 500 MS106, Batavia, IL 60510 



I. INTRODUCTION: EXPERIMENTAL STATUS AND SSM IMPROVEMENTS 



Since solar neutrinos were first detected two decades ago, the observed neutrino flux 
has always been a factor of 1.5 to 4 times less than that predicted by the Standard 
Solar Model (SSM). When both experimental and SSM uncertainties are included, the 
Homestake chlorine (CI) experiment, Kamiokande, and GALLEX rates respectively are 
approximately 6a, 3a, and 2a below the SSM predictions. Although the most recent 
GALLEX deficit is only 35%, significant differences from the SSM predictions persist. In 
this paper, we consider both astrophysical solutions (a cooler Sun model) and particle- 
physics solutions (MSW) to the solar neutrino problem. We show that no reasonable 
change in the SSM can reconcile the quoted Homestake, Kamiokande, and GALLEX 
results. On the other hand the MSW effect, which assumes neutrino mass differences 
Am 2 and mixing sin 2 29, accommodates all data: it does not require discarding any 
of the experiments, nor stretching the SSM beyond its uncertainties. Once MSW is 
admitted, the data determine T c at the 5% level, yielding a value consistent with the 
SSM. We also discuss possible MSW parameters when a non-standard core temperature 
is assumed. 

There are four measurements of the solar neutrinos so far. The chlorine experiment 
at Homestake [IJ is mainly sensitive to 8 B and 7 Be neutrinos, and reports^] an observed 
rate 2.1±0.3 SNU (0.28±0.04 of the SSM), while the SSM prediction is 8.0±1.0 SNU 
0. (Quoted errors are all 1 a in this paper.) A direct counting measurement by the 
Kamiokande Collaboration observes Cerenkov light from recoil electrons scattered by 
8 B neutrinos. The combined result of Kamiokande II J| and 220 days of Kamiokande III 
0] is 0.49±0.08 of the central value of the SSM; there is an additional 14% uncertainty 
in the SSM prediction. A unique opportunity to observe low energy pp neutrinos, 
which come from the main reaction responsible for the energy generation in the Sun, 
is provided by the SAGE and GALLEX gallium experiments. The GALLEX result is 
83±21 SNU (0.63±0.16 of the SSM) §, again well below the SSM value \32t\ SNU. 
SAGE has published a lower value 20±38 SNU based on their first 5 extractions ||. 
The later SAGE runs clearly suggest a higher rate, but no cumulative rate has been 
published. We will therefore exclude the SAGE results from our analysis. 

The SSM itself has undergone several refinements over the years. There are presently 

* We ignore the possibility of time dependence in the Homestake data. 
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at least four SSMs H 0, H §|, which, for the same physics input, agree with each other 
within 1% H and agree with the sound speed calculated in p-mode helioseismology 



within 2% QlO | . The latest published solar model by Bahcall and Pinsonneault improves 
on earlier calculations by using the most recent OPAL calculated opacities, meteoritic 
iron abundances, and updated nuclear-reaction cross sections M. The new model in- 
cludes the effects of helium diffusion so that the observed 3 He surface abundance is 
obtained. Their prediction for the convective zone boundary is in striking agreement 
with helioseismology data [TT]. The surface abundances of the 7 Li and 7 Be are still over- 
estimated, but the surface abundances have negligible effect on the neutrino production 
in the core. 



II. COOLER SUN WILL NOT EXPLAIN THE OBSERVED NEUTRINO FLUXES 

The production of high energy 8 B neutrinos and intermediate energy 7 Be neutrinos 
is directly proportional to the 3 He(a,7) 7 Be and 7 Be(p,7) 8 B nuclear cross sections and 
depends very sensitively on the solar temperature in the innermost 5% of the Sun's 
radius. This inner core temperature is not probed by existing p-wave helioseismological 
observations, but is determined by the radiative opacities throughout the Sun, which 
are believed to be calculated to within a few percent 0. 

Many nuclear and astrophysical explanations have been proposed to explain the 
solar neutrino deficit. One possibility is to change input parameters, such as lowering 
the 7 Be(p,7) 8 B cross section or reducing the opacity. Another is to invoke mechanisms 
that are not included in the SSM, such as a large core-magnetic field, a rapidly rotating 
core, or hypothetical WIMPs which carry away energy from the core. The net effect is 
generally a lowering of the core temperature and thus a reduction of the nuclear burning 
taking place there. 

It is therefore reasonable to examine the core temperature as a diagnostic of a whole 
class of astrophysical explanations of the neutrino deficit, although the temperature pro- 
file is the result, not an input, of a solar model [T2|]. We choose the central temperature 



(T c ) as a phenomenological parameter representing different conceivable solar models. 
The approximate correlation of the neutrino fluxes with T c was obtained by Bahcall 
and Ulrich |1J by examining 1000 self-consistent SSMs with randomly distributed input 
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parameters, and is given as simple power laws 0: 

0(pp) ~ r- u , 0( 7 Be) ~ T 8 , 0( 8 B ) ~ T 6 18 . (1) 

For each experiment we therefore parameterize rates relative to the SSM as functions 
of T c : 

R C i = (1 ±0.033) [0.775 x (1 ± 0.100) xT c 18 + 0.150 x (1± 0.036) xT c 8 
+ small terms] 
R Kam = (1 ± 0.100) xT c 18 
R Ga = (1 ±0.04) [0.538 x (1 ± 0.002) xT; L2 + 0.271 x (1 ± 0.036) x T 8 
+0.105 x (1 ± 0.100) x T c 18 + small terms], 

where T c is the central temperature relative to the SSM (T c = 1 = 15.67 x 10 6 K). The 
Bahcall-Pinsonneault solar model including diffusion M is used throughout the paper, 
unless otherwise mentioned. The uncertainty in the overall factors for CI and Ga is 
due to the detector reaction cross sections. The uncertainties in each neutrino flux 
include only nuclear physics uncertainties in production cross sections; astrophysical 
uncertainties are absorbed into a variable T c . These flux uncertainties are properly 
correlated for the three experiments. 

The T c dependence of the Kamiokande, CI, and Ga detectors are each shown in 
Fig. 1, and the best fits for various combinations of the data are summarized in Table I. 
If each experiment were fit alone, Kamiokande and Homestake require a reduction of 
T c by 4±1% and 10±1% respectively. For GALLEX alone, T c must be reduced by 
14±4%, but still does not completely fit the data, because the negative exponent of 
T c in the pp flux works against the reduction of the total rate. In order to match the 
observed luminosity, the pp chain has to compensate the energy loss due to reductions 
of the other reactions. As a result the T c reduction fails to fit the central value of the 
GALLEX rate, although it is compatible with the upper end of the range. The three 
separate T c fits are respectively 3a, 6a, and 3a below the SSM prediction T c = 1 ± AT C , 
where the SSM uncertainty AT C = 0.0057 is estimated from the 1000 SSM calculations 
in Ref. @. 

The combined observations cannot be fit by any single T c . The larger Kamiokande 
rate relative to the CI rate especially contradicts the T c dependence shown in Fig. 1. 



4 



The simultaneous fit of Kamiokande and CI yields T c = 0.92±0.01, but the x 2 value is so 
large (x 2 =13.78) as to exclude the fit at the >99.99% C.L. The combined Kamiokande 
and GALLEX results yield a marginally consistent T c : the best fit is T c = 0.96±0.01 
with x 2 — 2.64 (89.3% C.L.). When all three experiments are fit simultaneously, T c = 
0.92±0.01 but the x 2 test rejects the cooler Sun hypothesis at 99.99% C.L. This strong 
rejection of the cooler Sun is driven mainly by the contradiction between Kamiokande 
and Homestake and, secondarily, the low GALLEX result. 

In our fit, T c is allowed to vary in a range wider than that for which the power laws 
(Eqn. (1)) were derived. We assume that, even for large changes of T c , the power laws 
yield a reasonable approximation to the core-temperature dependence of the neutrino 
fluxes. It should be noted that our conclusions do not depend on a specific choice of the 
exponents given in Eqn. (1). Provided that the 7 Be flux is less temperature dependent 
than the 8 B flux, the CI rate is expected to be larger than the Kamiokande rate, 
contradicting the data [14"|. Even if both flux components had the same temperature 



exponent (=18), we find that T c = 0.93 ± 0.01 and x 2 = 10.21 for 1 degree of freedom: 
the x 2 test excludes the fit at 99.91% C.L., with similar conclusions for other exponents. 
Based on these observations one can conclude that the cooler Sun model, which is a 
nearly universal feature of astrophysical solutions to the solar neutrino problem, is 
strongly disfavored by the data: the central value of the GALLEX result and the 
combined result of Homestake and Kamiokande are each incompatible with the cooler 
Sun hypothesis. 

III. MSW FIT TO THE COMBINED OBSERVATIONS 

While modifications of the solar model cannot accommodate the data, an attractive 
solution is proposed from particle physics. Matter-enhanced neutrino oscillation (MSW 
effect), first proposed by Wolfenstein, then applied to solar neutrinos by Mikheyev and 
Smirnov, offers a natural explanation of the observed solar neutrino deficit without 
requiring any ad- hoc mechanisms |T5|]. This MSW mechanism assumes new properties 



of neutrinos, mass and mixings, to convert electron-neutrinos to other species when 
the neutrinos propagate through the Sun, making possible a large reduction of the 
^-counting rate. The MSW assumption of neutrino mass and mixing is a natural ex- 
tension of the Standard Model, and requires no ad-hoc features such as a large magnetic 
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moment. Unlike vacuum oscillations, it does not require fine-tuning. If the MSW os- 
cillation takes place in the Sun, the determination of the neutrino mass and mixing 
will provide a clue to grand unified theories, which naturally lead to parameters in the 
relevant region |TJ], [IT]]. 

The MSW effect depends on two parameters: the mass-squared difference Am 2 
and the vacuum mixing angle 9 between v e and another neutrino species into which it 
converts. The conversion occurs in a wide parameter space that covers four orders of 
magnitude both in Am 2 and sin 2 29: a triangle-shaped region in the sin 2 29 vs. Am 2 /E 
plane, surrounded by Am 2 /E < 2 x KT 5 (eV 2 /MeV) and sin 2 29 ■ Am 2 /E > 1(T 9 
(eV 2 /MeV), where E is neutrino energy. Typical survival-probability contours are 
shown in the Am 2 vs. sin 2 29 plane (MSW diagram) after integrations over the neutrino 



production site and the neutrino energy, including the detector cross sections |16|. 

The MSW diagrams show three physically-distinct regions: the adiabatic region 
(the horizontal upper arm of the triangle), the non-adiabatic region (the diagonal arm), 
and the large-mixing region (the right, vertical arm). For the adiabatic solution, the 
MSW resonance takes place in the core of the Sun where the neutrinos are produced; 
the density is high enough for the higher-energy neutrinos to resonate and be depleted 
while the lower-energy ones survive. For the non-adiabatic solution, the higher-energy 
neutrinos survive more because of non-adiabatic (Landau-Zener) jumping. In the large- 
mixing region, which connects smoothly to vacuum oscillation, the neutrino spectrum 
is equally reduced over the whole spectrum. In the middle of the isosnu triangle, almost 
100% conversion of v e occurs. 

This flexibility makes the MSW effect phenomenologically robust. It can preferen- 
tially suppress the high energy ( 8 B) neutrinos, or the low energy ( 7 Be and pp) neutrinos 
more. It can deplete the lower energy part of the 8 B and 7 Be spectrum while keeping 
the pp flux, as suggested by the experiments (Fig. 14). 

In the MSW calculations, we use the Parke formula [|1J|| instead of solving the 
Schrodinger equation for the oscillations numerically. The jump probability given by the 
exponential Landau-Zener formula Pj = e~ x is used [Til], where x = fthsm 2 9Am 2 /E, 
E is the neutrino energy, and h = (—d log n e /dr)~ 1 is the electron-density scale height 
evaluated at the resonance point. This formula agrees with the exact solution for large 



mixing angles and with the linear Landau-Zener approximation |20| for the small mixing 
region. The spatial distribution of neutrino production in the Sun is taken from Bahcall 
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and Pinsonneault || . The detector cross sections are taken from Bahcall and Ulrich . 
For Kamiokande, the (or u T ) contribution for flavor oscillations, the energy threshold, 
the energy resolution, and the trigger efficiency are all properly included |2T . 



In fitting the data, theoretical uncertainties of the SSM are treated with care, using 
a simple and transparent parameterization. For each flux component the nuclear cross 
section uncertainties of every reaction are added quadratically to the detector cross 
section uncertainties and to the astrophysical (non-nuclear) uncertainties. The latter 
is represented by the uncertainty in the central temperature AT C times the exponent 
defined in Eq (1). The theoretical uncertainty AT C = 0.0057 is chosen to yield flux un- 
certainties consistent with those given in Ref. 0, and to be consistent with estimates 
from the 1000 SSM Monte-Carlo calculations of Bahcall and Ulrich The corre- 

lations of the uncertainties among the experiments and flux components are properly 
taken into account. Our calculations were compared with other studies which utilize 



1000 Monte-Carlo SSMs |22|, |23f; the agreement is excellent. 

Effects of the astrophysical uncertainties on the MSW effect were also examined. 
Both the uncertainties from the neutrino production profile, which affects the matter 
mixing angle at the neutrino production, and from the electron-density scale height, 
which enters in the jump probability, were found to be small. 

Survival-probability contours and 90%-C.L. allowed regions are shown for each ex- 
periment in Fig. 2, 3(a), 4 and 5. The Homestake allowed region (Fig. 2) does not 
precisely trace the iso-probability contour because of the difference in 7 Be and 8 B 
theoretical uncertainties. Fig. 3(a) and 3 (b) show the Kamiokande result for flavor 
oscillations and oscillations to sterile neutrinos, respectively. 

The combined result (90% C.L.) of Homestake, Kamiokande II+III, and GALLEX 
is displayed in Fig. 6 (a). The confidence-level region is defined from x 2 values that sat- 
isfy x 2 (sin 2 29, Am 2 ) = Xmin + ^-Q, which is valid in the approximation that the allowed 
regions are 'ellipses' on the log sin 2 29 — log Am 2 plane. Including the GALLEX observa- 
tions, the allowed MSW parameters are either Am 2 = (0.3 — 1.2) x 10 -5 eV 2 , sin 2 29 = 
(0.4 - 1.5) x 10" 2 (non-adiabatic solution), or Am 2 = (0.3 - 4) x 10~ 5 eV 2 , sin 2 26* = 
0.5 — 0.9 (large-mixing solution). The best fits of Am 2 and sin 2 29 along with \ 2 value 
for each region are listed in the second and third columns of Table II. The experiments 
prefer the non-adiabatic solution to the large-angle solution. The non-adiabatic MSW 
solution yields a good fit (x 2 =0.56); but in the large-mixing region, the \ 2 value is 
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large (x 2 = 3.52): it is allowed at the 90% CL. by the definition above, but for 1 
degree of freedom (= 3 experiments — 2 parameters) this region is excluded at 95% 
CL. The allowed regions at the 68, 90, and 95% CL. are shown in Fig. 6(b): there 
is no parameter allowed in the large-angle region at 68% CL. The combined fit with- 
out theoretical uncertainties from the SSM and detector cross sections is displayed in 
Fig. 6 (c). Comparison with Fig. 6 (a) shows the noticeable effect of SSM uncertainties. 
Our Fig. 6(a) practically agrees with that obtained by the GALLEX group ||, who 
included the day-night effect and v e regeneration in the Earth, which we have neglected. 

Allowed regions for various combinations of any two experiments are shown in Fig. 7 
(Kamiokande and Homestake), Fig. 8 (Kamiokande and GALLEX), and Fig. 9 (Home- 
stake and GALLEX). 



We have also examined possibilities of oscillations to a sterile neutrino p4| , |25 
If v e oscillates into a sterile neutrino instead of or z/ T , the term n e — n n /2 enters 
the MSW equation in place of n e . (Here n e and n n are the local electron and neutron 
densities in the Sun ||24|| .) Also, for the Kamiokande detector there is no neutral current 
contribution from the converted neutrinos. The result for v e oscillations into sterile 
neutrinos is displayed in Fig. 3 (b) for Kamiokande, and in Fig. 10 and Table III for 
the combined fit. (There is no significant change for the Homestake and GALLEX 
experiments.) Because of the constraints by the Kamiokande observations, there is no 
solution in the large-angle region at 90% CL., and the allowed parameters are limited in 
the non-adiabatic region. Even in the non-adiabatic region, the best fit yields x 2 — 3.17 
and is excluded at 92% CL. for 1 degree of freedom. This is because of the smaller 
Homestake rate relative to Kamiokande: the absence of the neutral current events in 
Kamiokande requires a larger electron survival probability than for the flavor-oscillation 
case, and therefore widens the discrepancy between the two experiments. 

The precise determination of MSW parameters will allow us to draw some theoretical 
conclusions in Sec. V and to make predictions for next-generation neutrino experiments 
in Sec. VI. 

IV. SIMULTANEOUS FIT OF MSW AND CORE TEMPERATURE 

We have also studied the possibilities of having both MSW oscillations and a non- 
standard solar model by allowing T c to be a completely free parameter. (We use the 
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same nuclear and detector uncertainties as in the SSM case.) The data are fit simulta- 
neously to three parameters Am 2 , sin 2 29, and T c . The \ 2 pl°t is displayed as a function 
of T c (Fig. 11) , where x 2 is minimized for each T c with respect to Am 2 and sin 2 29 in 
the allowed region in the MSW diagram. By the y 2 fit, the data determine the core 
temperature at the 5% level. The best fits are T c = 1.03lo!o5 (90% C.L.) in the non- 
adiabatic region and T c = 1.05lo!o7 m the large-mixing region, in good agreement with 
the SSM prediction T c — 1 ±0.0057. The consistency between the data and the SSM is 
encouraging. Moreover, even allowing the MSW conversion and the other uncertainties, 
the observations determine the core temperature to within 5% [p6] . 

The allowed region for the 3-parameter fit is shown in Fig. 12. For each Am 2 and 
sin 2 29, the x 2 is minimized with respect to T c , and x 2 (sin 2 2$, Am 2 ) = Xmin + 4-6, 
determines the 90% C.L. allowed region, where Xmin i s a minimum with respect to all 
3 parameters. By allowing T c to be a free parameter, the allowed region is widened, 
now stretching from sin 2 29 = 4 x 10~ 4 ~ 1 and Am 2 = (0.2 — 4) x 10~ 5 eV 2 . The 
best fit parameters are shown in Tables II and III. The T c dependence of the region is 
seen in Fig. 13 (a) and Fig. 13 (b), which shows the 90% C.L. contours when T c is fixed 
at 1.05 and 0.95 respectively. The higher temperature Sun allows a region between 
the two islands allowed in the SSM case (Fig. 6(a)), while the cooler Sun pushes the 
parameter-space outward. 



V. THEORETICAL IMPLICATION OF FITTED MSW PARAMETERS 

The best fit MSW parameters from Fig. 6 (a) and Fig. 11 are summarized in Table II 
for the SSM (T c = 1) and for non-SSM in which T c is an adjustable parameter. The re- 
sults for sterile-neutrino oscillations are listed in Table III. The Am 2 range is consistent 
with the general expectations of grand unified theories WM or string-inspired models 



171 , but the mixing angles are not in agreement with the expectation 9i evtcm ~ 9qkm of 
the simplest GUTs |TJJ]. If we accept the SSM, the observed neutrino mixing are at least 
four times larger (or considerably smaller) than the quark mixings, sin 2 29 = 0.18 and 
< 2 x 10~ 3 for u — c and u—t quarks, respectively. If we allow a warm Sun, then Cabibbo 
mixing with Am 2 ~ (0.8 — 2) x 10 -5 eV 2 is possible; this suggests m v ~ (3 — 4) x 10 -3 
eV and an SO(10) GUT with intermediate-scale symmetry breaking. If we allow a cool 
Sun, then Am 2 ~ (0.5 — 1.5) x 1CT 5 eV 2 is possible, suggesting m Vr ~ (2 — 4) x 10" 3 eV 



9 



and a SUSY SO(10) GUT (with large Higgs representations so that the seesaw scale is 
close to the unification scale). 

GUT predictions for neutrino masses are much less robust than for mixing angles 
and mass ratios. We therefore regard the seesaw model [^7[] only as a crude guide to 
neutrino masses. In the SUSY GUT case, all neutrino masses are cosmologically and 
astrophysically insignificant. In the non-SUSY SO(10) GUT, the seesaw model suggests 
cosmologically interesting v T masses. Since m t /m c ~ 100 and m T /m M = 17 we have 
m„ T ~ 0.4 — 0.8 (40 — 80) eV for a linear (quadratic) seesaw model with up-quark masses 
and m„ r ~ 0.07 — 0.14 (1 — 2) eV for a linear (quadratic) seesaw model with charged 
lepton masses. In carrying out these extrapolations, we have taken m v — (2 — 4) x 10~ 3 
eV and, for the up-quark case, have included a factor two to renormalize mass ratio 
from from GUT to low energy scales [14]. (This factor increases beyond two non- 



linearly for large value of the top quark mass because of the Higgs corrections to the 
top mass PBfl.) We see that only a intermediate scale quadratic seesaw mechanism can 



give cosmologically significant v T mass [14]. String-inspired models can generate an 



intermediate seesaw scale via effective non-renormalizable operators [fTfl| . There are no 
clear mixing angle predictions in such models, but for v e — > oscillations consistent 
with either MSW solution, the v T may again be cosmologically significant. 



VI. PROSPECTS FOR FUTURE EXPERIMENTS 

One can predict the results of future solar neutrino observations from the parameters 
of the combined fit (Fig. 6 (a)). The v e survival probability is shown as a function of 
energy in Fig. 14 for each of the allowed regions. The predicted observed rates for the 
high-energy 8 B v e — e scattering (SNO and Super-Kamiokande 0]), the v e — d 
reaction (SNO), and the v e — e scattering from 7 Be neutrinos (BOREXINO [[H]]) are 
listed in Table IV |B3|. The detector cross sections as well as the proposed energy 



resolution of the detectors are included in the calculation. 

The measurement of the charged current reaction v e + d —> e + p + p, which is 
planned for the first-year operation of SNO, will be a clear diagnostic of the MSW; the 
distortion of the energy spectrum is characteristic of most particle physics solutions of 
the solar neutrino problem and cannot be caused by any astrophysical effects operative 
in the Sun ||33||. Figure 15 shows the predicted energy spectra for both the non-adiabatic 
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and the large-mixing solution, along with estimated statistical errors equivalent to a 
2- year operation (6000 total events). The non-adiabatic spectrum is very similar to 
the predicted spectrum for sterile neutrinos. If there is a non-adiabatic MSW effect, as 
suggested by the best fit (Table II), the spectral distortion will (a) confirm the MSW 
effect, and (b) discriminate between the two presently allowed solutions. Of course, the 
neutral to charged current ratio would also establish MSW oscillations into or u T 
(but not into a sterile v). 

The measurements of the neutral current events by SNO and BOREXINO, and 
their ratio relative to the charged current events would provide definite evidence for 
MSW flavor oscillations. The neutral current mode can determine the core temperature 
precisely even in the presence of the MSW effect: assuming 10% systematic errors, the 
high sensitivity of the 8 B neutrino flux (~ T c 18 ) allows a determination of T c at the 0.5% 
level. 

VII. SUMMARY 

Existing Homestake, Kamiokande, and GALLEX experiments strongly disfavor as- 
trophysical solutions invoking a cooler Sun. For the matter-enhanced neutrino oscil- 
lations, the data constrain the parameter space to two small regions, one in the non- 
adiabatic region (which is preferred by the data) and one in the large-mixing region. 
The fit for oscillations into a sterile neutrino allows only non-adiabatic oscillations at 
90% C.L. Allowing a non-standard core temperature along with MSW oscillations we 
find that the data constrain the core temperature at the 5% level, yielding values consis- 
tent with the SSM. The warmer Sun (T c =1.05) allows the parameter space predicted by 
the SO(10) GUT with an intermediate-breaking scale (for which the v T may be cosmo- 
logically relevant), while the cooler Sun stretches the allowed parameters into a region 
predicted by simple super symmetric SO(10) GUTs. Superstring-inspired models are 
consistent with all solutions. Predictions are made for future solar neutrino detectors. 
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model and neutrino fluxes in computer-readable form. It is a pleasure to thank Eugene 
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Table Captions 

Table I. The T c fits for various combinations of the Homestake (CI), Kamiokande II+III 
(Kam), and GALLEX results. Listed are the best fit value of T c with 90% confidence 
level (C.L.) error, x 2 values, and C.L. of excluding the fits. 

Table II. The best fit of the Homestake, Kamiokande II+III, and GALLEX results. In 
the T c — 1 (SSM) column, shown are the best fit of Am 2 and sin 2 29 with the x 2 value 
for each of the allowed MSW regions. The T c = free column is for the three parameter 
fit (Am 2 , sin 2 29, and T c ). The T c errors are at 90% C.L. 

Table III. The best fit of the Homestake, Kamiokande II+III, and GALLEX results for 
sterile-neutrino oscillations. In the T c — 1 (SSM) column, shown are the best fit of 
Am 2 and sin 2 28 with the x 2 value for each of the allowed MSW regions. The T c = free 
column is for the three parameter fit (Am 2 , sin 2 29, and T c ). The T c errors are at 90% 
C.L. 

Table IV. Predicted rates for future solar neutrino detectors, relative to the SSM expec- 
tations. The rates are listed for each of the allowed regions obtained from the best fit 
(Fig. 6(a)). 
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Figure Captions 

Figure 1. The approximate T c dependence of the neutrino counting rates (relative to 
the SSM) for the CI, Ga, and Kamiokande experiments, according to the power laws 
(Eqn. (1)). T c is relative to the SSM value (T c = 1 = 15.67 x 10 6 K = 1.35 keV). 

Figure 2. The v e survival probability contours (solid lines) for CI experiments and the 
allowed region obtained from the Homestake result (90% C.L., shaded region). The 
contours are for survival probabilities of 0.1, 0.2,- • •, 0.9, starting with the innermost 
solid line. The calculation of the allowed region includes the experimental errors, the 
detector cross section uncertainties, and the SSM flux errors. The allowed region slightly 
deviates from the iso-probability contours because of the difference in the 7 Be and 8 B 
flux uncertainties. 

Figure 3 (a). The v e survival probability contours (solid lines) for Kamiokande experiments 
and the allowed region for the Kamiokande II and III (220 days) result (90% C.L., 
shaded region). This is for flavor oscillations into or v T . The contours are for 
effective survival probabilities of 0.2, 0.3,- • •, 0.9, which include the effects of neutral 
current scattering. (There is no 0.1 contour.) The calculation includes the energy 
threshold, the energy resolution, and the trigger efficiency The allowed region includes 
the SSM uncertainties of the 8 B flux. 

Figure 3 (b). Same as Fig. 3 (a), except that it is for oscillations into a sterile neutrino. 
The contours are for survival probabilities of 0.1, 0.2,- • •, 0.9. Compared to flavor 
oscillations, lack of a neutral current contribution increases the v e survival probability 
required by the data and therefore pushes the allowed region outward of the triangle. 

Figure 4. The v e survival probability contours (solid lines) for Ga experiments and the 
90%-C.L. allowed region for the GALLEX result. The calculation of the allowed region 
includes the experimental errors, the detector cross section uncertainties, and the SSM 
flux errors. The contours are for survival probabilities of 0.1, 0.2,- • •, 0.9. 

Figure 5. Same as Fig. 4, except the allowed region is for SAGE. 

Figure 6 (a). The allowed region for Homestake (dotted region), Kamiokande II+III (solid 
line), and GALLEX (dashed line) results. The shaded region is the combined fit of the 
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three experiments (90% C.L.). 

Figure 6 (b). The allowed region of combined Homestake, Kamiokande II+III, and GAL- 
LEX results at 68% (shaded region), 90% (dashed lines), and 95% (dotted lines) C.L. 

Figure 6(c). The allowed region of the combined results with (dotted line) and without 
(shaded regions) theoretical uncertainties from the SSM and the detector cross sections. 

Figure 7. The combined fit of the Kamiokande II+III and Homestake results (90% C.L.). 

Figure 8. The combined fit of the Kamiokande II+III and GALLEX results (90% C.L.). 

Figure 9. The combined fit of the Homestake and GALLEX results (90% C.L.). 

Figure 10. The fit to oscillations to sterile neutrinos at 68% (shaded), 90% (dashed 
line) and 95% (dotted line) C.L. The data are the combined result of the Homestake, 
Kamiokande II+III, and GALLEX. 

Figure 11. The x 2 plot as a function of T c in the three parameter fit (T c , Am 2 , and sin 2 29) 
of the Homestake, Kamiokande, and GALLEX. For each T c , the x 2 is minimized with 
respect to Am 2 and sin 2 29. The data determine T c within 0.98 < T c < 1.06 (90% 
C.L.), which is consistent with the SSM value. 

Figure 12. The allowed MSW region of the combined result of Homestake, Kamiokande 
II+III, and GALLEX, using T c as a free parameter. As a result of allowing T c to change, 
the 90% C.L. region widens. Also shown are Am 2 , sin 2 29 predicted by SUSY SO(10) 
GUTs (shaded), intermediate-scale SO(10) GUTs (thick line), and string-inspired SUSY 
models with nonrenormalizable operators (shaded line). In each model the predictions 
for Am 2 are not robust and easy to change. In the string-inspired model, sin 2 29 is also 
changeable. 

Figure 13 (a). The allowed MSW region of the combined result of Homestake, Kamiokande 
II+III, and GALLEX when T c is fixed at 1.05 (a warmer Sun). As a result of the high 
T c , the two regions of the SSM fit (Fig. 6 (a)) merge into one, allowing the parameters 
predicted by SO(10) GUTs with intermediate-breaking scales. 

Figure 13 (b). The allowed MSW region of the combined result of Homestake, Kamiokande 
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II+III, and GALLEX when T c is fixed at 0.95 (a cooler Sun). As a result of the low 
T c , the allowed region is pushed outward of the MSW region, allowing the parameters 
predicted by supersymmetric SO(10) GUTs. 

Figure 14 The u e survival probabilities as a function of energy for the two regions obtained 
by the best fit (Fig. 6 (a)). The solid line is for the non-adiabatic region and the dashed 
line is for the large-mixing region. 

Figure 15. The predicted spectral shape of charged-current events at SNO. The solid 
line is for the allowed parameter space in the non-adiabatic region, and the dashed line 
is for the large-mixing region. The errors are equivalent to a 2-year operation (6000 
events). The distortion of the spectrum in the non-adiabatic branch will confirm the 
MSW effect and differentiate the two allowed regions (Fig. 6 (a)). 
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Table I: T c Fit for Kamiokande II+III, Homestake, and GALLEX 





T c ± AT C 


Y 2 
A 


G L 

V . J_J . 


Kam 


0.961 


± 


0.010 







CI 


0.901 


± 


0.015 







GALLEX 


0.860 


± 


0.042 


1.31 




Kam+Cl 


0.921 


± 


0.013 


13.78 


>99.99 


Kam+GALLEX 


0.960 


± 


0.010 


2.64 


89.26 


Cl+GALLEX 


0.902 


± 


0.013 


1.49 


77.56 


Kam+Cl+GALLEX 


0.920 


± 


0.012 


15.46 


99.99 



Table II: Best Fit of Kamiokande II+III, Homestake, and GALLEX 





T 

± c 

Non-adiabatic 


= 1 
Large- mixing 


T = 

Non-adiabatic 


^free 
Large-mixing 


sin 2 29 

Am 2 (eV 2 ) 

T c 

2 

X 


6.7xl0~ 3 
6.1xl0- 6 
1+0.0057 
0.56 


0.73 
l.lxl0~ 5 
1+0.0057 

3.52 


7.7xl0~ 3 
8.9xl0~ 6 
1 03+ - 03 

i - uo -0.05 

0. 


0.31 
l.lxl0~ 5 

i - ud -0.07 

0. 



Table HI: Best Fit for Oscillations into Sterile Neutrino 





T 

± c 

Non-adiabatic 


= 1 
Large- mixing 


T = 

Non-adiabatic 


;free 
Large-mixing 


sin 2 20 

Am 2 (eV 2 ) 

T c 

x 2 


9.7xl0~ 3 
3.8xl0~ 6 
1+0.0057 
3.17 


0.83 
4.8xl0~ 6 
1+0.0057 

8.95 


7.5 xl0~ 3 

4.0xl0- 6 
n qq+ - 16 

u - yy -0.02 

3.13 


0.15 

3.5xl0~ 6 

i i cr+0.005 
1 - io -0.18 

4.91 
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Table IV: Predicted Rates for Future Detectors 



SNO (Charged Current) 
SNO (y — e scattering) 
Super-Kamiokande 
Borexino ( 7 Bez/ — e) 



Rate / SSM 
Non-adiabatic Large-mixing 
0.15-0.5 0.15-0.3 
0.3 - 0.6 0.3 - 0.45 

0.3 - 0.6 0.3 - 0.45 

0.15-0.6 0.3-0.6 
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